Data Analysis of Continuous Gravitational Wave Signal: Fourier Transform 
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We present the Fourier Transform of continuous gravitational wave for arbitrary location of detector 
and source and for any duration of observation time in which both rotational motion of earth about 
its spin axis and orbital motion around sun has been taken into account. We also give the method 
to account the spin down of continuous gravitational wave. 
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The detection of gravitational waves (GW) from their 
possible sources, at the present stage, has to largely de- 
pend upon the careful study of the expected characteris- 
tic of potential sources and the waveforms. In this con- 
text, the continuous sources, in particular pulsars, are 
now receiving much attention [1,2]. These may be the 
one of the sources to look for the signals in data streams 
of the long arm laser interferometer and are also attrac- 
tive sources for currently operating bar detectors. The 
theoretical research clearly shows that the data analysis 
will involve processing of very long time series and ex- 
ploration of a very large parameter space. It turns out 
that optimal processing of month long data streams is 
computationally prohibitive [1] and methods to reduce 
the computation burden are being worked out [4,5]. 

At the root cause of the excessive computational re- 
quirement lies the fourier transform (FT) analysis of the 
data. In order to detect the signal from the dominant 
noise one has to analyse the long time observation data 
ranging from months to year. The output of detector will 
be Doppler modulated; both frequency and amplitude 
due to the rotations of earth. The maximum Doppler 
shift for one kHz signal is ~ 21.42 x IQ-^ [8]. The out- 
put depends in a complex manner on the location of the 
source and the detector. Consequently, the FT is ob- 
tained using numerical techniques; one commonly in use 
is fast fourier transform (FFT). The evalutaion of FT an- 
alytically is an important step in the data analysis and 
in contrast to FFT has following advantages: 

• The computational time is extremely less and par- 
ticularly in the case of long data train gains signif- 
icance. 

• We can achieve the resolution of the transform as 
we need, hence dominant peaks of doppler modu- 
lated signals can be accumulated, and 

• To make templates for matched filtering in colored 



In the literature very little attention has been paid to 
analytical analysis for fourier transform [6-8] . In this let- 
ter we address to this problem in its complete generality. 



The response R{t) of the ground based laser interfer- 
ometric detector at time i is a linear combination of the 
two polarizations ft.+ and ft.x of the signal given via 



R{t) ^ F+{t)h+{t) + Fy,{t)hy,{t) 



(1) 



where F^ , F-^ are the beam pattern functions and repre- 
sent the amplitude modulation (AM) part. In Solar Sys- 
tem Barycentre (SSB) frame the beam pattern as func- 
tion of the direction of the incoming wave (9,(j),ip), the 
orientation of the detector angles (a, /?, 7) and e, the an- 
gle between the equator and the ecliptic may be expressed 
as [6] 

F+ = i [2{L^ - M^)ZU + (N^ - P^)A - 

(Q2 - R'^)B] + {LN - MP)C + [LQ + MR)D 



F^ 



where 



+ {NQ + PR)E, 

2LMZU + NPA - ]-B sin^ 9 sin 2(1) 

{LP + MN)C + [MQ - LR)D 
+ {PQ - NR)E 



(2) 



(3) 



A = 2XY cos^ e — sin^ e sin^ a sin 27 -|- 

sin2e(Xsina sin 7 — Y sin a cos 7) , 
B = 2XY sin^ e — cos^ e sin^ a sin 27 — 

sin 2e{X sin a sin 'y — Y sin a cos 7) , 
C — cos e{YV + XU) + sin e(y sin a sin 7 — C/ sin a cos 7) , 
D = — sine(yy + XU) + cose(y sina sin7 — [/sin a cos 7] , 
E — —2XY cos e sin e — cos e sin e sin^ a sin 27 -|- 

cos 2e{X sin a sin 'y — Y sin a cos 7) , 
L = cos ip cos (f> — cosO sin <j) sin ?/; , 
M = cos V' sin (j) + cosO cos (j) sin ij: , 
N = — sin -0 cos (j) — cosO sin cos ■0 , 
P = — sin ip sin -|- cosO cos cos ip , 
Q — sin 9 sin , R = sin 9 cos 4> , 

[/ = — cos a cos (3 sin 7 — sin (3 cos 7 , 
V — cos a cos P cos 7 — sin /3 sin 7 , 



X = COS a sin (3 cos 7 + cos f3 sin 7 , 
Y = — cos a sin /3 sin 7 + cos /3 cos 7 , 

(3 — Po + Wrt 

Here /3o represent the initial azimuth of the detector. We 
have assumed that the arms of the laser interferometer 
are mutually orthogonal. However, if the angle between 
the arms of the detector are different than 90° the expres- 
sion for F^ and F^ will simply get multiplied by sine of 
the angle between the arms of the detector respectively. 
Inclusion of this parameter does not pose any difficulty 
as it simply affect the norm of the FT by the same factor. 

After algebraic manipulation eq. (2) and (3) may be 
expressed as 

F+ = Fi^ cos 2/3 + F2+ sin 2/3 + F^^^ cos/3 + 

f4+sin/3 + F5^ (4) 

F^ = Fi.^ COS2/3 + F2, sin 2/3 + Fg^ cos/3 + 

F4, sin/3 + F5, (5) 

Fi^ and Fi^ [i ~ 1,2,3,4,5) are time independent and 
are given by 



Fi, = —20 cos a cos 27 



iJ sin 27 , o , 

(cos a + 1) . 



2 

F2^ = H cos a cos 27 + G sin 27(cos^ a + 1) , 
F^^ = I sin a cos 27 + J sin 2a sin 27 , 

I 

Fi^ ~ 2 J sin a cos 27 — — sin 2a sin 27 , 

3 sin a sin 27 



F5 



-[H + L^ - M^], 



1 



G = - [{LQ + MR) sin e - [LN - MP) cos e] 

H = ]: [{N'^ - P^) cos^ e-{L^ - M^) + 

(g2 - R^) sin^ e - (NQ + PR) sin 2e] , 
I ^- [(Q2 _ ^2^ gjj^ 2e - {N^ - P^) sin 2e - 

2(7Vg + Pi?)cos2e], 
J = - [{LN - MP) sin e + {LQ + Mi?) cos e] 



and Fi^ is related to Fi^ via 



-Fjx (^^ 'P^ V', a, /?, 7, e) = ^j+ (^, </> - ^, ^, «, /?, 7, e) (6) 

i = 1,2,3,4,5 

The two polarisation states of the signal may be taken as 

h+{t) = ho^cos[^{t)], h^{t) = ho^sm[^{t)] 

where ft-o+j ^ox ^^'^ the time independent amplitude of 
h+{t), and hx{t). Here $(t) is phase of the GW signal 
and may be expressed in SSB frame as [7] 



$(t) = 27rfo 



R. 



Rs 



■ sm W cos ( 



— sin aisin 0(sin /3 cos e sin </> + cos (f) cos /3) 
c 



sin /3 sine cos 0} 
R. 



R. 



■ sin 9 cos ( 



sin ajsin 0(sin /3o cos e sin + cos (j) cos /5o) + 
c 

sin /3o sine cos 0}] (7) 

where, </>' — Worbt~4'i Rse is the distance from the centre 
of the SSB frame to the centre of the earth, i?e is the ra- 
dius of the earth, Wgrh is the angular velocity of the earth 
around the sun, Wr is the angular velocity of the earth 
about its spin axis and c is the velocity of light. Here Jo 
is the frequency of the continuous signal. 

The complete response of the detector is given via 

R{f) = R+{I) + RAJ) , (8) 

R+{f) = jh+{t)F+{t)e-^^^f'dt, (9) 

i?x (/) = /" /ix {t)F^ {t)e-'^^f'dt (10) 

Due to the symmetries involved in F^ and F-^ as given 
by [eqn. (6)] it is sufRcient to evaluate either of the FTs 
(9) and (10) and the other may be obtained in a simple 
manner. We will, at present, consider R^{f) and get 



R+if) 



0*2/3, 



-i2/3c 



(Fi++*i^2+)/i+(./ + 2/,)/2+ 



e*^'^°(Fl^-^F2J/l+(/-2^)/2 + 

{F3^+zFi^)h+{f + fr)/2 + 



'tPa 



^iPc 



(i^3+ - iFi^)h+{I - ./r)/2 + F^^h{f) (11) 



where 



h+{I) 



ho+ cos[$(t)] 



(12) 



and represent the FT arising due to frequency modula- 
tion. This means that the AM results in the four side 
bands at frequencies / ± /^ , / ± 2/^ around the central 
frequency /. For the moment, we keep the data interval 
of the observation time completely arbitrary. We evalu- 
ate the FT 

An+l)At 

h{f) = h+{.f) = / Vcos[<I'(t)] (13) 

JnAt 

where n is an integer and At is time interval. It is im- 
portant to point out that it is not easy to evaluate the 
FT analytically and it is the stage where one has to take 
recourse to numerical methods. The FT [eq. (13)] after 
lengthy but straight forward calculation is obtained as a 
double series 



kf) 



2Wr 



fc— oo ?n— oo 



Y^ J2 e^M2M)[C{XX,r) 



-cxD m— — oo 



iV{X,C,t)] 



(14) 



where 



^ - (/o - /)//,. , 
Jfe(Z)J„(A/-) 



(/r = Wr /2Tr) , 



kX 



<) 



C(A, C, r) = sin(z^T) cos(afcr + ttit — fcA — to(^) 

ak + m N • / , 
|cos(z/TJ sin(aKT + rriT — 

V 

+ sin(fcA + toC)} 7 

2?(A, C, t) = cos(z^r) cos(afcr + rtiT — fcA — to^) 

afc + m 
H sin(i/rj siii(afcr + ?7ir — kA — ?7i{^j 

V 

— cos(fcA + toC) , 
7^ =: Z cos + Q + 27r/on A t , 

X = (j) — anr , ( = 6 — nr , 

S — tan^ — ■ , T 

Q 
V — (27r/c)/oi?e sin a [cos /?o(sin 9 cos e sin + 

cos sin e) — sin [3o sin cos (j))] , 
Q = {2tt / c) fo Re sin a [sin /3o (sin 6* cose sin;/) + 

cos 6 sin e) + cos /3o sin 9 cos 0] , 

Jk{Z), Jmi-M) represent the bessel function of first kind 
of integer order k and m respectively. 



w^At, a ==1/365.25, 



The result obtained [eq. (14)] is useful in getting the 
FT of the complete response of a detector of GW be- 
ing emitted by a pulsar whose frequency is either slow- 
ing down or picking up. To account for this aspect, we 
consider the evaluation in different window of time by 
splitting the interval (0, Tobs) in N equal parts, each of 
interval At [Tabs — N At) such that the signal over a 
window may be treated as monochromatic. The strategy 
is to evaluate the FT over the window and finally to add 
the result. This process has also been suggested by Brady 
& Creighton [6] and Schutz [4] in numerical computing 
and called by them as stacking. The transform [eq. (14)] 
can be computed easily for data of any observation time. 
Consequently, the obtained transform may be one of the 
efficient tool for the detection of continuous gravitational 
wave. To obtain the result one has to evaluate R{f) in 
each window by accounting for changing value of /„. 
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A case of particular interest arises where one needs to 
evaluate the FT of the data for Tabs time i.e. the limits 
in (14) are taken as to Tobs- The result in this case is 
obtained by taking 



n = , At = Tobs 

Thus we obtain 
ho. V 



WrT 



r-^ obs 



Kf) 



2'Wr 



J2 E e'M2,J^)[C{XX,WrTob 

k= — oo m = — oo 
iV{X,C,WrTobs)] 



(15) 



Suppose one is interested in evaluating the FT for a day 
or year observation data then he/she can obtain the FT 
by putting n = in eq. (14), and 



Tobs = one year 
Tobs = one day 



= 27r/a 
27r 
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